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Abstract 

A body of literature has developed concerning "cloaking by anomalous local- 
ized resonance". The mathematical heart of the matter involves the behavior of 
a divergence-form elliptic equation in the plane, V • (a(x)Vtt(x)) = f(x). The 
complex-valued coefficient has a matrix-shell-core geometry, with real part equal 
to 1 in the matrix and the core, and -1 in the shell; one is interested in understand- 
ing the resonant behavior of the solution as the imaginary part of a{x) decreases 
to zero (so that ellipticity is lost). Most analytical work in this area has relied 
on separation of variables, and has therefore been restricted to radial geometries. 
We introduce a new approach based on a pair of dual variational principles, and 
apply it to some non-radial examples. In our examples, as in the radial setting, 
the spatial location of the source / plays a crucial role in determining whether or 
not resonance occurs. 

MSC: 35Q60, 35P05 

Keywords: cloaking, anomalous localized resonance, negative index metama- 
terials 

1 Introduction 

A body of literature has developed concerning "cloaking by anomalous localized reso- 
nance" . Cloaking of two types of objects has been considered: (a) dipoles or inclusions, 
considered e.g. in [2,4,8,10,13] and (b) spatially localized sources [1]. In this article we 
work in the setting of localized sources. 
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¥ u ' B q (0) 



B R (0) 

Figure 1: Sketch of the core- shell-matrix geometry. 



Focusing initially on the math (not the physics), we are interested in a divergence- 
form PDE in the plane: 

V-{a v Vu v ) = f onR 2 (1.1) 
Vu v — > as \x\ — > oo. 

The coefficient a v (x) is piecewise constant and complex- valued, with constant imaginary 
part r] > 0: 

a v (x) = A(x) + in; (1.2) 
its real part has a matrix-shell-core character in the sense that 

{+1 outside B R (0) 
-1 in the shell B R (0) \ E (1.3) 
+1 in the core S 

(see Figure 1). Concerning the core, we assume that 

£ C £i(0) (1.4) 

so the shell includes an annulus of width R — 1. Concerning the source /, we assume it 
is real-valued, supported at distance q from the origin, and has zero mean: 

f = FV}\dB q {Q), F : dB q (0) ->■ R, F e L 2 (dB q (0)), and / F alU 1 = . (1.5) 

JdB q {0) 

Our interest lies in the question: 

Question: yls — )• wit/i / and A(x) held fixed, what is the behavior of 

E v-=l i \Vu v \ 2 dxl 

^ ./1B2 
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In the radial setting, i.e. when E = -Bi(O), one expects by analogy with [1,2,8,10,13] 
that the answer depends mainly on the location of the source. Specifically: there is a 
critical radius R* = R 3 ^ 2 such that for a broad class of sources /, 



Note that it is no restriction to fix the core radius to be 1. A scaling argument implies 
that, for core radius r , and shell radius R, the critical radius is 



Definition 1.1 (Resonance). Let a configuration be given by coefficients A and source 
f as in (1.3)— (1.5) . We shall call the configuration resonant , if 



Otherwise we call the configuration non-resonant . 

In the physics literature the term "anomalous localized resonance" is used. An 
anomalous feature of the resonance is that it is not associated to a finite dimensional 
eigenvalue of a linear operator and a forcing term at or near the resonant frequency. 
Instead, the resonance here is associated to an infinite dimensional kernel of the limiting 
(non-elliptic) operator. The word localized refers to the fact that the resonance is spa- 
tially localized: while J |Viir,| 2 — > oo if q < R*, the potential u v (and therefore also its 
gradient Vu v ) stay uniformly bounded outside some ball. 

The connection to cloaking is as follows (see [1] for a more thorough discussion). 
For time-harmonic wave propagation in the quasistatic regime, E n is the rate at which 
energy is dissipated to heat. Let us now consider a source a v f, where a v G K is a scaling 
factor. If the (unsealed) source / produces resonance (i.e. if E v — > oo) then the source 
a v f is connected to the energy dissipation a^E^. If the physical source a v f has finite 
power, then we must have a v — > as r\ — > 0. If the fields u v associated with the unsealed 
source / are bounded outside a certain region, then the physical fields a. n u v vanish in 
that region as r\ — > 0. This implies that the finite power source a n f is not visible from 
outside. 

Returning to a more mathematical perspective: we are interested in this problem 
because it involves the behavior of the elliptic system (1.1) in a limit when ellipticity is 
lost. It is not surprising that oscillatory behavior occurs in such a limit. It is however 
surprising that, at least in the radial examples, (a) resonance depends so strongly on the 
location of the source, and (b) the oscillatory behavior is spatially localized. We would 
like to understand the following question: 

Question: Is this surprising behavior particular to the radial setting, or is it a more 
general phenomenon"? 



limsupi^ = oo if q < R*, while 



limsup E v < oo if q > R*. 




lim sup E v = oo . 
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The present paper addresses only point (a): the dependence on the location of the 
source. Our method, which is variational in character, is unfortunately not well-suited 
to the study of point (b). 

We know only one numerical study of a similar problem with non-radial (and non- 
slab) geometry. The paper [4] by Bruno and Lintner considered, via numerical simu- 
lation, various examples including an elliptical core in an elliptical shell. The results 
were similar to those of the radial case; in particular, the structure seemed to cloak a 
polarizable dipole placed sufficiently near the shell. 

The paper [1] by Ammari et al considers a problem very similar to ours. The main 
difference is that both the outer and inner edges of the shell are not constrained to 
be radial. (There is also a minor difference: their PDE has a v = 1 in the matrix and 
core and a v = — 1 + in in the shell, so energy is dissipated only in the shell.) Using a 
representation based on single layer potentials, Ammari et al obtain an expression for 
a spatially localized analog of E v . To make use of their expression, one needs detailed 
information on the spectral properties of certain boundary integral operators. This 
information is difficult to come by in general and hence, beyond the radial setting, it is 
unclear how to use their method to obtain information on resonance and non-resonance 
in the limit as rj — > 0. 

Our approach is based on variational principles. The starting point is a pair of (dual) 
variational principles for E v . One expresses E v as a minimum; trial functions may be 
used to provide an upper bound in order to show that resonance doesn't occur. The 
dual principle expresses E v as a maximum; trial functions may be used to provide a 
lower bound in order to show that resonance occurs. Similar variational principles were 
considered in [5,11]. Our main results - all proved using the variational principles - are 
the following: 

(i) If there is no core then there is always resonance, for any source radius q > R and 
any nonzero / (see Proposition 3.2). 

(ii) For any core E C -Bi(O), there is resonance for a broad class of sources /, provided 
the source location is q < R* := R 3 ^ 2 (see Theorem 3.4). 

(iii) In the radial case (when E = £?i(0)), R* := i? 3 / 2 is critical, in the sense that (a) 
when the source location is q < R* resonance occurs for a broad class of /'s, and 
(b) when the location is q > R* resonance does not occur for any / (see Theorem 
3.4 and Proposition 4.1). 

(iv) In the (weakly) nonradial case when the core is B p (zq) with |z | sufficiently near 
and p sufficiently near 1, resonance does not occur if the source location q is 
sufficiently large (see Theorem 5.3). 

Point (iii) is already known, from Section 5 of [1] . Our variational method is interest- 
ing even in this radial setting: our proof of (iii) is, we think, simpler and more elementary 
than the argument of [1]. Unfortunately, our methods do not seem to provide simple 
proofs for the localization effect when there is resonance. 

In focusing on (1.1), we have chosen the imaginary part of a v to be the same con- 
stant constant in the matrix, shell, and core. This simplifies the formulas, and it seems 
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physically unobjectionable. But we suppose a similar method could be used when the 
imaginary part is different in each region. 

We specifically consider sources / that are concentrated on the curve dB q (0). Our 
method also allows the study of more general distributions of sources, which can be 
obtained as a superpositions of concentrated sources, /, at different values of q. 

We have taken the core to have A(x) = 1 because this case has particular interest: 
in the radial setting, the "cloaking device is invisible" if the core has A — 1, see [13]. 
However anomalous localized resonance also occurs when A takes a different (constant) 
value in the core. It would be interesting to extend our method to analyze cores with 



Our assumption that A{x) = —1 in the shell is essential to the phenomenon. Indeed, 
our PDE problem becomes very different if the ratio across each interface, the plasmonic 
eigenvalue, is different from —1. This can be seen from the perspective of the boundary 
integral method, where ratios other than —1 lead to boundary integral equations of 
Fredholm type, see [1,7] 

Our main results are almost exclusively for a circular outer shell boundary dB R (0). 
This is essential to our method, since we use the perfect plasmon waves on the outer 
shell boundary in the construction of comparison functions. We refer to Section 3 for a 
further discussion. A more general geometry is only treated in Proposition 3.3 with the 
help of a domain transformation. Related techniques are used in [12]. 

Plasmonic resonance effects have many potential applications. This is one of the 
reasons why the development of negative index metamaterials is another much-studied 
research area, see e.g. [3,9,14]. We hope that our variational approach will be useful 
also in the other contexts. 

Notation. We use polar coordinates and write x G M? as x = r (cos 6, sin 9). In 
Section 5 we identify M 2 = C via (x±, x 2 ) = X\ +ix 2 = z. With this notation, we identify 
z = re ld . The complex conjugate of z is denoted by z. 

We denote the sphere with radius p, centered at x , as B p (x ). The measure T-^ydVt 
is the 1-dimensional Hausdorff measure on the curve dVt. Unless otherwise specified, 
integrals are over all of M 2 . Constants C may change from one line to the next. 

2 The primal and dual variational principles 

In the subsequent definitions of energies we always consider the source / as a given 
element / G H^ 1 {W'). We will always consider sources with a compact support (in the 
sense of distributions). Furthermore, we shall assume that the sources / have a vanishing 
average, 



Since / is merely a distribution, it would be more correct to write (/, 1) = 0, where 
1 : R 2 — > R is the constant function, l(x) = 1 for all x G M 2 . We note that since / has 
compact support, it can be applied to test-functions that are only locally of class H 1 . 



A^l. 
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We remark that, while the main results of this paper concern R 2 , the primal and dual 
variational principles generalize to any dimension. 

2.1 A complex elliptic system and its non-elliptic limit 

Our aim is to study, for a sequence 77 = 77^ — >■ 0, sequences u v of solutions to (1.1). 
For non- vanishing dissipation, 77 7^ 0, (1.1) is an elliptic PDE, while the system loses 
ellipticity in the limit 77—7-0. 

To a solution u v : R 2 — > C of the original complex-valued equation 

V-(a„V<> = / (2-1) 

we have associated an energy E v (in physical terms the energy dissipation in the struc- 
ture) 

*,K):=? / |V«,| 2 . (2.2) 

As noted in the introduction, the phenomenon of cloaking is related to resonance in the 
sense of Definition 1.1, 

E v (u v ) 00 (2.3) 

along a subsequence 77 \ 0. 

We can write the complex scalar equation for u v : R 2 — > C as a system of two real 
scalar equations. We set 

u v = v v + i — w v , with v v , w v : R 2 — > R. (2.4) 

For a real-valued source, / : R 2 — > R, the complex equation V • (a v Vu v ) = f with 
a v = A + in is equivalent to the coupled system of two real equations on R 2 , 

V • (AVv v ) - Aw v = /, (2.5) 
V • {AVw v ) + r] 2 Av v = 0. (2.6) 

The energy E v (u v ) can be expressed in terms of v v and w v as 

* f 1^.. |2_ V f ,2 , 1 f ,V7„„ 12 



^K) = \ j |V«,| 2 = ^ J \Vv v \ 2 + -J \Vw n \ 2 . (2.7) 
In the following subsections we introduce 

1. the primal variational problem, a minimization problem, which characterizes the 
energy E v {u^) as a constrained minimum; and 

2. the dual variational problem, a maximization problem, which characterizes the 
energy E v (u v ) as a constrained maximum. 

To provide a functional analytic framework for the study of the variational problems 
we introduce the following function space of real or complex-valued functions, 

^(R 2 ) := {U E L 2 oc (R 2 ) \VU E L 2 (R 2 ) } , \\U\\% := [ \VU\ 2 + [ \U\ 2 . 

(2.8) 
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2.2 The primal variational problem 

For fixed / G if _1 (IR 2 ) we consider the energy functional 



defined for v,w G if (R ). The primal variational problem is given by 

minimize I v (v,w) over all pairs (v,w) 

which satisfy the PDE constraint V • (AVv) — Aw = f . 



(2.10) 



Lemma 2.1. Let f G ff^IR 2 ) &e a fixed real-valued source with compact support and 
with vanishing average. Then the primal variational problem (2.10) is equivalent to the 
original problem (2.1) with energy (2.2) in the following sense. 

1. The infimum 

inf | I v (v,w) (v,w) G H 1 (M 2 ) x H 1 (M 2 ), V • (AVv) - Aw = /} (2.11) 

is attained at a pair (v v ,w v ) G ii^M 2 ) x if^M 2 ). 

2. The minimizing pair, (v v ,w v ), is unique up to an additive constant. The function 
u v := v^ + irj^Wrj is the unique (up to an additive constant) solution of the original 
problem (2.1). 

3. For the solutions, the energies coincide, 

E v (u v ) = iw„) . (2.12) 

Remark. The lemma implies 

E v (u v ) < I v (v,w) (2.13) 

for every pair (v,w) that satisfies the PDE constraint of (2.10). We shall use the in- 
equality (2.13) to establish non-resonance. 

Proof. Point 1. Fix a radius s > such that supp(/) C B s (0), we introduce the function 
space with constraint: 

X := lueHH 



u = 

B s (0) 

Note that I v , defined in (2.9), is convex on X x X. Moreover, the constraint set is non- 
empty. Indeed, choose v* smooth and of compact support and defined w* to be the weak 
solution of Aw t = — / + V ■ (AVv m ). It follows that the infimum in (2.11) is attained on 
X x X (see, e.g., Chapter 8.2 of [6]), i.e. there exists (v v ,w v ) G X x X such that 

In(v v , w v ) < I v (v, w) for all (v, w) G X x X, with V • (AVv) — Aw = / . 
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Point 2. We first observe that the constraint (2.10) is identical to (2.5). As a minimizer 
of Ijj, the pair (v v ,w v ) satisfies the Euler-Lagrange equation 

d T I v ( v v + tv, w v + tw ) = 0, for every (v, w) G X x X 

V / r=0 

satisfying V • (AVv) — Aw = 0. For the energy I v , this equation reads 

= 7] J Vv v ■ Vv + i J Vw v ■ Vw = 7] J Vv v ■ Vv + i J Vw r , ■ AVv, 
where we have used the constraint to obtain the second equality. We find 
-^(rfA Vrt + V ■ {AVw v ),vj = 0, Vtiel, 

which is the weak form of (2.6). We use here that v can be any element of if 1 (R 2 ) 
with compact support, since an associated w G ij 1 (R 2 ) can be obtained as the solution 
of a Poisson problem. As a solution of (2.5)-(2.6), the pair (v v ,w v ) defines through 
u v = v r) + vr\~ x w n a solution of the original problem (2.1). 

The uniqueness is a consequence of the fact that the original problem (2.1) possesses 
a unique solution. This can be seen from the Lax-Milgram Lemma. We introduce a 
sesquilinear form b(-, •) : X x X — > C defined by 



u 2 • 



b(u\,U2) := — i / a v VuiV 
The form b(-, •) is coercive on X 

m(u,u) > C\\u\\ 2 x , 

since the imaginary part of a v is strictly positive and by the Poincare inequality. Ex- 
istence and uniqueness of a weak solution u v G H 1 ^ 2 ) solution follows from the Lax- 
Milgram Lemma. 

Point 3. The energy equality (2.12) was already observed in (2.7). 

The proof of Lemma 2.1 is complete. □ 

2.3 The dual variational problem 

For fixed / G if" 1 (M 2 ), we introduce the dual energy 

JM) ■= |/V - I J - I J |W| 2 , (2.14) 

defined for (v,ip) G if^M 2 ). The dual variational problem is given by 

maximize J v (v,ip) over all pairs (v,ip) 
which satisfy the PDE constraint V • (AVtp) + nAv = 0. 



(2.15) 



The following lemma establishes that the dual variational problem is also equivalent 
to the original complex equation. 
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Lemma 2.2. Let f G iJ _1 (]R 2 ) &e a /ixec? real-valued source with compact support and 
with vanishing average. Then the dual variational problem (2.15) is equivalent to the 
original problem (2.1) with energy (2.2) in the following sense. 

1. The supremum 

supj J v (v,i>) (v,i>) G H^M 2 ) x H\R 2 ), V • (AVi>) + i]Av = oj (2.16) 

is attained at a pair (v v ,ip v ) G if^R 2 ) x i7 1 (M 2 ). 

2. The maximizing pair, {v v ,ip v ), is unique up to an additive constant. The function 
u v := Vrf + iipr! is the unique (up to constants) solution of the original problem (2.1). 

3. For the solutions, the energies coincide, 

E v (u v ) = Jr,(Vr,,1pri) ■ (2-17) 

Remark. The lemma implies 

> J v {v,$) (2.18) 

for every pair (v, ifj), which satisfies the PDE constraint of (2.15). We shall use inequality 
(2.18) to establish our results on resonance. 

Proof. Point 1. The existence of a maximizing pair (i^,^) f° r t ne variational problem 
(2.16) follows from the concavity of J v (convexity of — J v ), by arguments analogous to 
those given above for the primal variational problem. 

Point 2. At a maximizer, (y v , ip v ), one has for all v : ip G if 1 (lR 2 ) fl L 2 (R 2 ) satisfying the 
PDE constraint V • (AV$) + r]Av = 0, that 

d T J v (vri + tv, ip v + ripj = 0. 

For the energy J n , this relation provides 

J ftp-r] J Vv v -Vv-r] J VVvW* = 0. (2.19) 

Using the PDE constraint V • (AVtp) + r/Av = to replace v, we find that (2.19) is 
equivalent to 



= J ' flp + J VV V ■ AVlp -7} J Vtpr, ■ Vip 

= (/ - V • (AVu„) + tjA^) . 



We conclude that the pair (v v ,w v ) := (^,77^) is a weak solution of (2.5)-(2.6) and, 
thus, that u v := v v + iip v is a solution of V • (a^Wujj) — f on M. 2 . Uniqueness follows 
again from the fact that u v is unique up to constants. 
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Point 3. Regarding the energy equality, we calculate 

E^Ur,) - J^Vr,,^) = | J | VwJ 2 - J„(u„ ^) = V J \Vv v \ 2 + 1] J \Vl^ v \ 2 - J f^ v 
= - J Vu„ ■ AVlP, - (ifir,, -f + V • (AVVr,)) - J flp v = 0. 

This concludes the proof of Lemma 2.2. □ 

3 Resonance results 

As discussed in the introduction, we consider configurations of the following type: 

1. The coefficients a v (x) and A(x) are defined by (1.2)-(1.3) with core £ C £>i(0) 

2. The source, f(x), is concentrated at a distance q > from the origin and is taken 
of the form / = FH l [dB q (0) as in (1.5). 

We seek conditions on configurations, which ensure resonance or non-resonance in the 
sense of Definition 1.1. 

We explore the resonance properties of a configuration as follows. To prove resonance 
we use the dual variational principle, exploiting (2.18). It suffices to construct, given 
rj = rjj — > 0, a sequence of comparison functions {v v ^ v ) that satisfy the constraint of 
(2.15) and that have unbounded energies J^v^,^). To prove non-resonance we use the 
primal variational principle, exploiting (2.13). It suffices to construct, given i] = rjj — >■ 0, 
a sequence of comparison functions (v v ,w v ) that satisfy the constraint of (2.10) which 
have bounded energies I^ivj^w^). 

In this section, we show resonance in both radial and non-radial settings. The non- 
resonance results will be presented in Section 4 for radial cases and Section 5 for a 
non-radial geometry. 

The basis of construction of trial functions is the family of perfect plasmon waves: 
Remark 3.1. Consider the PDE for functions ip : IR 2 — > R, 

V-(AV^) = 0, 

Vij}{x) ->■ 0, as \x\ ->■ oo (3.1) 

where 

m = {~) (3-2) 



For any k > 1, there is a non-trivial solution ip = ipk which achieves its maximum at a 
point with \x\ = R, given by: 



M*)--={ rk ?t 6) f° r \ x \< R > (3.3) 

1 r .-k R 2k cos ( fc 0) f or \ x \ > Rm 
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We call such functions perfect plasmon waves. Notice that 

[ |V^.| 2 = / ij k 

J J\x\=R 

Since our proofs rely on these perfect plasmon waves, we always use (except for 
Proposition 3.3) the circular outer shell boundary dB R (0). For the same reason, our 
technique is restricted to the two-dimensional setting. (For an explanation why it does 
not extend to 3 or more dimensions, see Appendix A). 

Since perfect plasmon waves are given in terms of Fourier harmonics, it is natu- 
ral to expand arbitrary sources F E L 2 (dB q (0)) with vanishing average in Fourier se- 
ries. We will represent an arbitrary source by superposition of the elementary sources 
parametrized by harmonic index, k EN, and source- distance, q: 

fl = cos{kff)'H 1 \dB q {{S). (3.4) 



dr 



2nkR 2k 



3.1 Resonance in the radial case 

Proposition 3.2 (No core =>- Resonance for sources at any distance q). Assume no 
core, £ = 0, so that a v (x) = A(x)+ir) where A(x) is given by (3.2). Let f = FH 1 [dB q (0) 
with 7^ F : dB q (0) — > K. be a source at a distance q > R. Then the configuration is 
resonant, i.e. E v (u v ) — >■ oo as rj — > 0. 

Proof. We fix the radii R and q and consider an arbitrary sequence rj — rjj — > 0. We 
write the source as / = YlT=i a kfk^ where f\ is defined in (3.4). Since F ^ 0, there 
exists some k > 1, such that a k ^ 0. Our aim is to find a sequence (w^,^), satisfying 
the constraint V • (AVip v ) + i]^ v v = of (2.15) and such that J ri {y v ,ip r ^ — > oo. We 
choose 

v v (x) = (3.5) 
ip v (x) := X v ipk(x) , (3.6) 

where is the perfect plasmon wave of (3.3) and A r; E K. is to be chosen below. The pair 
(v n ,ip v ) satisfies the constraint (2.15). Using (2.18), the definition of J v , the hypothesis 
q > R, and the orthogonality of Fourier harmonics, we obtain 

E v (u v ) > J ri (v v , 4> v ) = J v (0, ij) v )= J f ■ ip v - | J | V^| 2 

= / a k cos{k6) ■ X n q- k R 2k cos{k6) -^W 2 I |V^ fc | 2 

JdB q (0) 2 J 

> 7iqa k X v q~ k R 2k -C (r)\\ v \ 2 ) kR 2k . 
Choosing A^ — > oo with 77 1 | 2 — > we obtain E v (u v ) — >■ 00 for 77—7- 0. □ 
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3.2 Resonance in the non-radial case 

Our first observation for non-radial geometry regards a variant of Proposition 3.2. We 
consider the index A = —1 in a domain DcK 2 , which is similar, but not identical to 
the ball B R (0). 

Let $ : C — > C be a holomorphic function. For three radii R < q < s we introduce 
the three domains D R := $(B R (0)), D q := $(£,(())) and D s : = $(5 S (0)). We assume 
that $ is bijective on the largest ball, §\b s (o) '■ B s (0) — > D s has a inverse \& : D s — )■ B s (0). 

Proposition 3.3 (Resonance for non-radial structures without core). Fix radii 1 < R < 
q < s, holomorphic maps $ and \& = $~ 1 as above. Assume s > q 2 /R. Consider the 
equation (1.1).' 

V • (a^Vwr,) = / 

where a v (x) = A(x) + in and 

Then there exists a source f = Fl-i l \dD q , where F e L 2 (dD q ), such that the configura- 
tion is resonant, i.e. E v {u v ) — > oo for r\ — > 0. 

Proof. We proceed as in Proposition 3.2 and exploit the dual variational principle. Our 
aim is to construct a sequence (v v ,ip v ), satisfying the constraint V • {AWip v ) +i]Av v = 
of (2.15) and such that J v (v v ,ip v ) — > oo. However in this case, due to the coupling of 
Fourier harmonics in a non-radial geometry, we cannot restrict ourselves to a harmonic 
of fixed index, k. 

We start the construction from the perfect plasmon waves ipk of (3.3). They are 
mapped with $ to functions ipk '■ D s — > R, ^ := ^ o We note that these functions 
are harmonic in D s \ dD R . Regarding the jump of the normal derivative along dD R , 
we can calculate as follows. The normal vector e r to dB R (0) is mapped to the (not 
normalized) normal vector v = ■ e r of dD R . In a point x E dD R the matrix _D$ 
can be represented by a complex number M G C, and we can calculate (u, Vipk) = 
(Me r , ((V^ fc ) T • M~ l ) T ) = (Me r , (M" 1 ) T V^ fc ) = (e r , V^ fc ). Thus, ^ k solves 

V • (A(x)Vip k (xfj = 0, for x E D s . (3.8) 

In order to define functions on all of M 2 , we introduce a smooth cut-off function 
h : M 2 — > [0, 1]. By our assumption on the radii, we can choose a number Q E (q 2 /R, s) 
and D Q = $(B Q (0)). Let h be such that h = 1 on D Q and h = on M 2 \ D s . 

We now construct a comparison function of the form 

ipr,(x) := \ v h(x) tp k (x) , (3.9) 

where and k = k(r/) are to be chosen below. Once ip v is determined, the function v v 
is chosen as the bounded solution of 



V Av v (x) = -V • (A(x)Vi> v (x)) = -V • (A(x)V^(x)) 1 d ,\ Dq (x) . (3.10) 
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The latter equality holds by (3.8) since Vh(x) ^ precisely on D s \Dq. The pair (v v , ip v ) 
satisfies, by construction, the constraint (2.15). It only remains to verify Jn{v v , ^r?) ~~ ^ °°- 
To motivate our choice of A^ and k{rj), we first compute the contributions to the 
energy, J v , of the functions v v and ip v . In the following calculations, C denotes a constant 
that is independent of k and r\\ C depends on the geometry of the structure and its value 
may change from one line to the next. For the contribution from v v we find 



V J 



Vv n \ z <c- 



<c 



V Jd s \d q 



V Jd 3 \d q 

The contribution to the energy, J v , from ip v is 



\H 2 + Wk\ 2 <c^k 

V 



XI fR 2 ^ 2k 



Q 



(3.11) 



7/ 



J iv^l 2 < 



Cr]X 2 I kR 2k + 



1 



2\ 2fc N 



k \ Q 



2 r>2k 



< Cr]kX 2 R 



(3.12) 



It remains to evaluate the first term of the energy, J v , for a source / = F^^dDg. For 
some positive constant c > we obtain 



I fip v = X v [ F^ k >cX v I (Fo $)^ fc 

J JdD q JdB q (0) 

= Mt)7 < F °* )cos(M) = cA ' (?) 



Oik , 



(3.13) 



where a k is a Fourier coefficient for F o $. We shall drive J^v^^ip^) to infinity by 
driving the contribution (3.13) to infinity while keeping the contributions (3.11) and 
(3.12) bounded. 

Balancing the upper bounds (3.11) and (3.12) requires that we choose k = k{rj) such 
that: 



(R/Q) 



k{rf] 



rj. 



Since R/Q < R 2 /q 2 < 1, there is such k(r)) : and k(rf) — > oo as rj | 0. With this choice of 

k, 

?3\ Hv) 



/{■ / d3 \ ' 

|V^| 2 + 77 J |V^,| 2 <CX 2 v k(r))l^j 

To keep this contribution to the energy bounded we choose so that 



1 



\2 . 

" k( V ) \R* 



Q 



Hi) 



(3.14) 



(3.15) 



Finally, substitution of (3.15) into (3.13) we obtain, for some c, C > 0: 



> 



Mr, 



C > c 



QR\ 



k/2 



) 



a k 



c. 



(3.16) 
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We have chosen Q with Q > q 2 /R. Therefore, if we choose F such that its Fourier 
coefficients, a k , are not decaying too rapidly, we have | J^v^, i^ v )\ — > oo as rj — > 0. This 
completes the proof of Proposition 3.3. □ 

Next, we consider a non-radial geometry with a core, E C -Bi(O), of arbitrary shape. 
The following resonance result has a proof very similar to that of Proposition 3.2. 

Theorem 3.4 (Any shape core, source at q < R* =>- Resonance). Let E C -Bi(O) 
be an arbitrary core with F = <9E a curve of class C 2 . Then, for every radius R < 
q < R* := R 3 / 2 , there exists f supported at a distance q, such that the configuration is 
resonant. 

Remark. Our proof actually gives a slightly stronger result. Not only does there exist 
/, supported at q G (R,R*), E v — > oo as rj — > 0, but furthermore the divergence of 
the energy occurs for every / having high frequency components of sufficiently large 
amplitude. 

Proof. We fix R < q < R* and a sequence rj — rjj — > 0. We consider the source function 
/ = Sfeli a kfk with fl as in (3.4). Our aim is to construct a sequence (v v , ip v ), satisfying 
V • (AVip v ) + r]Av v = of (2.15) with J v (v v ,^ v ) ->■ oo. 

As in the proof of Proposition 3.2 our sequence of trial functions is built using perfect 
plasmon waves; as before we choose the harmonic index k = k(rj) to depend on r\ and 
set 

The numbers k = k{rf) G N and A,, G R will be chosen below. 

The function ip v is not A-harmonic along the core interface <9E C -Bi(O). In order to 
satisfy the constraint we therefore define v v as the solution of rj Aiv, = —V ■ (AVip v )- By 
elliptic estimates 

»7l|VuJi 2(R2) < Cn- 1 ||V- (AV^)||i-i (M2) < Cf 1 X 2 v k( v ). 

It remains to calculate the energy J v (v v ,ip v ). We choose k = k(n) to be the smallest 
integer with R~ k < rj and note that R~ k+l > rj holds. Exploiting (2.18), we obtain, for 
some Co > 

e v (u v ) > j„k,^) = y m - \ J iw.i 2 - \j i v ^i 2 

> c a k \ v q~ k R 2k - Cnk\ 2 v R 2k - Crf^k^ 
= \ v R k ^ c a k (^j - C\ v k(riR k ) - C-^X.k 

The choice of k with 1 < r/R k < R ensures that the last two contributions are of 
comparable order. We find, for some C > 0, 

/ /R\ k< - v) 
E v (u v ) > \R k ^ \c Q a m (-) -C X v k(r]) 
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We choose such that the right hand side is positive, specifically 



1 fR\ k{v] 



Thus, 



2 CoakM UJ J = 4C^) {CoakM)2 {-J J 

By assumption, q, the location of the source satisfies q < R*, or equivalently g 2 < R 3 . 
To ensure that E v (u v ) — > oo it suffices to assume that the sequence of Fourier coefficients 
((Xk)k ^ / 2 (N) decays sufficiently slowly. In particular, if (otk)k decays algebraically we 
have E v (u v ) — > oo. This completes the proof of Theorem 3.4. □ 



4 Non-resonance in the radial case 

In this section, we use the primal variational principle (2.10) to show non- resonance in 
the radial case for sources located at distance larger than the critical radius R*. 

Proposition 4.1 (Non- resonance beyond R* in the radial case). Consider a n (x) of 
(1.2)— (1.3), with the radial concentric arrangement E = -Bi(O) C B R (0). Assume f = 
FU l [dB q {Q), F E L 2 (dB q (0)). Then, for any q > R* := R 3 / 2 , the configuration is 
non-resonant. 

Proof. Expanding in Fourier series we have F = J2k>i a k cos (k@) + J2k>i Pk sin(k9) = 

F even + -^odd- It Suffices to prOVC that / cvcn = Feven'H 1 \pB q (0) and /odd = F ddH l [dB q (0) 

are non- resonant. We give the argument for / ev en; the argument for / dd is the same. 
Accordingly, we consider from now on / = J2k a kf k ' wnere /fe is given by (3.4) and 
(c*fc)fc>i e / 2 (N;1R); we suppress here the superscript q of f%. We will construct test 
functions in Step 1 and compute their energy in Step 2 to prove the Proposition. 

Step 1. Construction of comparison functions. To prove non- resonance, we use the 
primal variational problem (2.10). Consider a fixed sequence r\ = rjj tending to zero. We 
shall construct (v v ,w^), satisfying the constraint 

V-(AVv n ) - Aw v = f (4.1) 

such that the energy along this sequence, In(v v , w v ) remains bounded. Our strategy is 
to decompose the source / into a low frequency part and a high frequency part as 

k* oo 

f = /' ow + / high , / low :=^a fc / fc , / hish := £ a k f k , (4.2) 

k=l k=k*+l 

where k* is chosen to depend on rj. Later we will choose k* = k*{rf) to be the smallest 
integer for which R~ k * > rj. 
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Our approach, to be discussed in detail below, is to solve the constraint equation 
(4.1) in the form: v v = v l ™ + v^ igh where 



v L ° w satisfies V • (AV^ ow ) = f ow (4.3) 
^satisfies V • (AV^ igh ) = / high (4.4) 

w v satisfies - Aw n = -V • (AVv^ igh ) + / high (4.5) 

This construction yields (v v ,w v ), which satisfies the constraint (4.1) of the primal prob- 
lem (2.10). Furthermore, we shall see that with an appropriate choice of cutoff k* = k*(n) 
in (4.2), I^v^jW^) remains bounded as rj — > 0. As in our analysis of resonance, we shall 
make strong use of the perfect plasmon waves. 

Step la. Construction of v l ° w . The function v l ° w is pieced together using variants of the 
perfect plasmon waves. 



v k (x) := < 



( r k cos(k9) for |x| < 1, 

r~ k cos(k9) for 1 < \x\ < R, 

r k R~ 2k cos(k9) for R < \x\ < q, 

j~ k (q/R) 2k cos(k8) for q < \x\. 



(4.6) 



We note that v k has the following properties. 

1. v k is continuous on all M 2 

2. v k satisfies V • {AVv k ) = for x e M 2 \ dB q (0). 

3. Along dB q (0), v k has a jump in its normal flux: 



W ■ Vv k ] dB ( o) = {—q k R~ 2k - -q k R- 2k \ cos(^) = -— q k R~ 2k cos(k0) 
{ q q J q 



Therefore, an appropriate constant multiple \ k v k will satisfy 

V • (AV\ k v k ) = a k f k onR 2 . 
In order to satisfy this relation, we must choose with 

2k 



\/,- \-jq k R- 2k ) 



We therefore set 



v ^ := J2\ k v k , with \ k :=-a k ^-q- k R 2k . (4.7) 



k=i 
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This function satisfies (4.3), 

V • (AW™) = [u • W™] dBq{0) H 1 l_dB q (0) = A* (-—q k R~ 2k cos{k9)\ U 1 [_dB q (0) 
= ^a fc cos(A;^ 1 L^(0) = / 1 ' 



fc<fc* 

flow 



Step lb. Construction of t^ lgh and w v . The function v^ lgh is constructed from the 

elementary plasmon waves V k for the radius q. The functions are not tuned to solve 
V • (AVv) = on dBi(0) or dB R (0), but they are small along these curves (compared 
to their maximal values). We set 



%{x) : = 



r k cos(k9) 
r~ k q 2k cos(k9) 



for \x\ < q, 
for q < \x\. 



(4.8) 



Recall A(x) = 1 in a neighborhood of \x\ = q, so the jump in the normal flux on dB q (0) 
is 



v ■ AW k 



dB q (0) 



Therefore if we set 



k>k* 



-2k/q)q k cos(k6) . 



x Q -k 

k := ~ ak 2k q 



(4.9) 



it follows that (4.4) is satisfied: 

v • (Av^; igh ) 



dB q (0) 



f 



high 



We emphasize that v ^ lgh is not a solution on all of R 2 due to normal flux jumps at 
and \x\ = R. Since (4.3) is satisfied, the constraint (4.1) is equivalent to (4.5), 

-Aw v = -V • (AV^ igh ) + / high 



x\ 



-J2 Xk 



v ■ AVV k 



k>k* 



afii(o) 



-H}\dB x (G)-Y,^ 



v ■ AVV k 



k>k* 



dB R (0) 



W}[dB R (Q) 



(4.10) 

We use this equation to define w v . 

Step 2. Calculation of energies. It remains to calculate the energy I n {y v ,w^), for the 
above choice of v v = v l ° w + t>^ lgh and w v . It is in this step that we choose the low-high 
frequency cutoff, k* = k*(r)) to ensure that I v (v v ,w v ) remains uniformly bounded as 
77 — >• 0. Once we verify the boundedness of this sequence of energies, the non-resonance 
property of Definition 1.1 follows from (2.13). 

Step 2a. Energy of v v . By the triangle inequality we can bound the energies of v 1 ™ and 
v^ lgh separately. Furthermore, orthogonality of Fourier modes implies for v l ° w : 

?2 x 2k 



v / iv<° w i 2 = v J2 i Afc i 2 / i v ^i 2 ^ Cr > \ ak 



k<k* 



k<k* 



q 



k 2 



max 



(4.11) 
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For the case where q > R 2 , we obtain 

V [ \Vv l ™\ 2 <C V J2\<Xk\ 2 <Cri : 

J k<k* 

which is obviously bounded. The case where R* < q < R 2 is more subtle. We note here 
that estimate (4.11) simplifies in this case to 

v [W^Cn^Hl 2 (^T • (4.12) 
J k<k* v y / 



We will come back to this bound soon with a specific choice of k*. 
The energy of t>^ lgh is easier to control: 

V [ I V< is T <C V J2 l^l 2 / |VT4| 2 <C V J2 Kl' < C. (4.13) 
J k>k* J k 

Step 2b. Energy ofw v . Next we study the energy of w v . By the properties of the solution 
operator (— A) -1 acting on functions in if _1 (R 2 ) with vanishing average, we have 

- f \Vw v \ 2 < C-\\V ■ (AVv^ h ) - f hish \\ 2 H-^ <C-Y \\ k \ 2 R 2k k. 

The last inequality follows from (4.10), which states that V-(AVf^ lgh )— / hlgh is supported 
on \x\ = 1 and \x\ = R. 

Now by the choice of X k in (4.9), we have \X k \ < C\a k \k~ 1 q~ k , and hence 

2k* 



- [ \Vw v \ 2 < C £ Kl 2 ~ (-) ( 4 - 14 ) 



Balancing the right hand sides of the bounds (4.12) and (4.14) we choose k* so that 

R 2 \ 2k * 1 {R^ 2k * 



q J rj \ q 

i.e. we choose k* = k*(r)) to be the smallest integer with R~ k * < r\ such that 

r] < R- k * +l , and - < R k \ (4.15) 

V 



Combining (4.15) with (4.14) and (4.12), we obtain 

i J iv^p < c y, w 2 (£) (4.i6) 



and 



nj\W™\ 2 < c£K| 2 (^J (4.17) 



Thus, if q > R* = R 3 ^ 2 , Iiv^w^ is bounded as r\ — > 0. The proof of non-resonance is 
complete. □ 
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5 A non-resonance result in a non-circular geometry 
5.1 Interaction coefficients 

In this section we use complex notation. We will use complex analysis in order to 
calculate certain interaction integrals that will be of interest for non-resonance results 
in non-radial geometries. 

We identify 1R 2 = C via (xi,^) = x i + 1x2 = z. The complex functions z ^ z k and 
z h-> z~ k are holomorphic on C and C \ {0}; we have used the real part of these functions 
before, 

^ z k ) = r k cos(Jfe0), ^(z~ k ) = r~ k cos{k6), for z = re ie . 

By the Cauchy-Riemann differential equations, the gradient of these real functions can 
easily be calculated, 



VK(z k ) = (£) U{z k ) = (-J (z k ) ee d Xl (z k ) = (z k )> = kz k -' . 

Accordingly, we can evaluate for v = z/\z\ the normal derivative on a sphere dB r (0), 

{v, V^(z k )) = §R ( • kz k ~ l \ = -^(z k ) = -r k cos(k9) . (5.1) 
\\z\ J r r 

A similar calculation provides for the imaginary part (i/,V3(2 fc )) = j^s(z k ). 

Using the complex notation we define, for arbitrary radius p > and arbitrary center 
z G C, the function 

vjy J K ((*-*o) m ) for I^-^qI <p, 

mW " \p 2 "^((^-^ )- m ) for |^-^o| >P- 1 } 

This function is a perfect plasmon wave for dB p (zo) in the sense that it is continuous and 
it satisfies the equation V • (A ZOiP V^ m ) = for the coefficient A ZQjP = 1 outside B p (z ) 
and A ZOjf) = — 1 inside B p (z ). This is easily checked with either the above complex 
calculations or the previously imployed real calculations. 

Of importance will be the interaction of two perfect plasmon waves with different 
centers. In particular, we need to know the following interaction coefficients. 

Definition 5.1. Let a radius p > and a center zq G C be such that the corresponding 
disk contains the origin, G £ := B p (z ). For wave numbers k,m GN, the interaction 
coefficient is defined as 

I m , k := I M(z- k )(z - z ) m dH 1 GC. (5.3) 

JdT, 

The integral in (5.3) is with respect to the Hausdorff measure; the complex number 
Im,k can therefore be identified with the real vector whose two components are the 
real integrals J $t(z~ k )$t(z — z ) m dH 1 and J $t(z~ k )?s(z — zq) 171 dV} . In this sense, the 
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coefficients I m ^ are, up to normalizing factors, the coefficients for an expansion of the 
function $l(z~ k ) in spherical harmonics of the sphere dB p (zo). 

With the help of complex analysis, the interaction coefficients can be computed ex- 
plicitly. 

Lemma 5.2 (Properties of I m ,k)- The interaction coefficients for m > 1 and k>0 are 

given by 

{0 for m < k, 

i i\m-k (m-l\ m „ fc ^ , (5.4) 

(-l) m k np\ ] z form>k 

\m — k J 

For m = we have io,fc — 27rp 5o,fc- 

On the circle dB p (z ), the function $t,(z~ k ) with k > 1 can be expanded in spherical 
harmonics with center z as 

^( Z - k ) = ^Y,^p 2m+1 y lI ^ k ■ ( z - z ^ m ■ ( 5 - 5 ) 

m>k 

For any number Q > we have the estimate 

Y^Q k \^A < vpQdzol + QT 1 - 1 (5.6) 

ken 

for every m G N. 

Proof. We note that, for m = k = 0, the value 7 ,o = 27rp follows immediately from the 
definition of l m ,k- From now on, we can therefore assume k + m > 1. 

We can calculate the number with the help of the residue theorem. We decom- 
pose the integral according to 3t,(z~ k ) = \z~ k + \z~ k . We calculate first the contribution 
of z~ k , 

( [ z- k (z - z ) m dn l ) = [ z- k {Y^) m dU 1 = p 2m [ z- k (z - z )- m dH 1 

® -ipp 2m f z~ k (z - z,y m ' 1 dz ( = } , 

JdY, 

where the symbol c.c. denotes complex conjugation. In the calculation above, we have 
used in equality (a) the fact that, for the argument i? G K. of (z — z ), we have 
(z — zo) m = p m exp(—im$) = p 2m (z — z )~ m . In equality (b) we have introduced the 
complex line element dz with the help of the tangential vector i(z — z )/p, substituting 
(z — zq) dH 1 = —ipdz. In equality (c) we have used the fact that the contour of integra- 
tion can be deformed without changing the value of the integral; we deform the contour 
into increasingly large circles, and the limiting value is 0. We exploited here k + m > 1. 

We have thus seen that one of the two contributions to I m ^ vanishes. Using again 
the tangential line element [z — z ) dH 1 = —ipdz we have 

Im,k = \ I Z~ k { Z - Z ) m dV} = ^ I Z~ k (Z - Z ) m ~ l dz . 
1 JdT. 1 JdT, 
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At this point we have verified the claims for m = 0, since for k > 1 we can again move 
the contour of integration to oo, hence the integral vanishes. In the case m > 1 we 
expand the term (z — Zq) 111 ^ 1 and find 

Im, k = £ ("J ') dz. 

By the residue theorem, the boundary integral is non-vanishing only for j = m — k, since 
for this value the exponent of z is ra — 1— j — k — — 1. We obtain 



- =r L z " t - i ) <-*>"-* * - ( : : I ) <-*> 



m—k 



by Cauchy's theorem. This proves the explicit formula (5.4). 

For fixed A; > 1, we expand the function $l(z~ k ) [dB p (z ) : dB p (z ) — >■ R in spherical 
harmonics with coefficients 7/, 7; G R, 

= ^( z - Z °Y + ^( z - z oY • (5-7) 

For arbitrary m > 1, we multiply this function with 3?(z — ^ ) m and integrate over 
T := dB p (z ). Using orthogonality properties of spherical harmonics, we find 

/ ^{z- k )^(z - z ) m dU 1 = J2li zo) l M(z - z ) m dH 1 = i m irp 2m+1 . 

The left hand side is nothing else than dt(I m% k)- Repeating the calculation for the imag- 
inary part, we find 7 m + i^ m = (irp 2m+1 )~ 1 I m) k. This verifies the expansion (5.5). 

Estimate (5.6) is obtained with a straightforward calculation. For m > 1 there holds 

Y,Q k \^A< E ^(Z-k) \ z o\ m -~ k Q k = ^pq (\ z o\ + Qr- 1 ■ 

kefi l<k<m ^ ' 

This completes the proof. □ 



5.2 Non-resonance for a non-concentric circular core 

The following result generalizes Proposition 4.1 to a geometry that is not radially sym- 
metric. 

Theorem 5.3 (Non-resonance for non-concentric core). We consider a configuration 
of the following form. The coefficients are given by a circular core X = B p (z ) with 
G £ C -Bi(O) through (1.2) -(1.3). The source is located at a radius q with q > R 3 , and 
given as f = Fl-i^ydBqifS) with F = ^2 k>1 ol^ cos(A;6>) + ^ fc>1 fik sm(k9). We assume 
that the Fourier coefficients satisfy ^ k {\(Xk\ + \Pk\} — C . 

There exist Eq = Eo(q) > and E\ = E\(q) > such that, if the core is close to the unit 
disk in the sense that \z \ < e and 1 — E\ < p < 1, the configuration is non-resonant. 
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Remark. We will provide an explicit condition regarding the smallness of 6$ and e±, see 
(5.14) and (5.19). 

Note that we consider sources at the radius q > R 3 . We know that the source radius 
q must satisfy q > R* = R 3 ^ 2 to be non-resonant, see Theorem 3.4. The lower bound R 3 
is probably not optimal. 

Proof. As in the proof of Proposition 4.1, we can decompose the expansion of F into two 
parts and can write F = Feven + ^odd- By linearity of the equations it is sufficient to show 
the non-resonance property for the two contributions separately. Without restriction of 
generality, we study in the following / = J2 k a k fk with f k = /| of (3.4), and coefficients 
(a k ) k G l 2 (N; R) fl / 1 (N; R). 

We fix a sequence r) — rjj \ 0. Our aim is to construct a sequence (v v , w v ) of bounded 
energy I v , and to use the primal variational principle (2.10) to show non-resonance. 

Step 1. Construction of comparison functions. We will use a construction similar in 
spirit to that used in the proof of Proposition 4.1. The main difference is that the 
functions v k (x) of (4.6) are not suited for the eccentric core E. We will replace these 
functions by v k (x) of (5.8) defined in Step la. We then need to correct errors on <9E due 
to the non-concentric geometry. This is done in Step lb and Step lc. 

Step la. Construction of the main part of v v . We first construct the main part of v v , 
denoted by V v , following the construction of Proposition 4.1 with the new elementary 
functions 

'v k (x)\x for i6E, 



v k (x) :-- 



r 



-k 



cos(k9) for x E B R (0) \ E, 



(5.8) 

r k R~ 2k cos(k9) for R < \x\ < q, 

r~ k {q/R) 2k cos(k8) for q < \x\, 

where Vk\i: is chosen so that Vk is harmonic in E and continuous on <9E. 

The function V v is constructed as a linear combination of the elementary functions 
Vk- The coefficients are chosen to satisfy V • (AVV^) = f away from 9E. This leads 
to 

V n := ^ X k v k , A, := -a k ^-q- k R 2k . (5.9) 

km 

The coefficients are actually identical to those in (4.7). This is because Vk coincides 
with v k on R 2 \ E. 

Step lb. Evaluation of errors on <9E. In the case of concentric spheres, the construction 
could be finished at this point, the distinction into high and low frequencies was only 
necessary in order to find the optimal bound for q. Instead, since we now study a core E 
that is not concentric, the functions v k are not solutions of V • (AVv) = on <9E. Hence, 
we need to correct the error on <9E: 



F := V • (AVV V ) -f = J2^ RvfcLt + cUfcU U 1 L<9E 

km 

= X kduVkU U l L«9S + AfcWin U l L^E (5.10) 



km keN 

F^H^d^ + F^yd^, 
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where the last equality gives the definition of F out and F m . 

We start with F out , the contributions from d u Vk\ ut, which can be calculated explicitly, 
since Vk is defined in (5.8) as Vk(z) = r~ h cos(k6) = 3l(z~ k ) on x G -Br(O) \ E. As in (5.1) 
we calculate the normal derivative, which we then expand using (5.5), 

m>k 

Hence, we have evaluated the first part of the error in (5.10) to be 

E Xk ^^2 » • (* - z °D = E R W • (* - ^) m ) ( 5 - n ) 



out 

^ — ' 7rp^ 

k,m>l m>l 



with 



Kfc<m ^ 



We next estimate the decay of asm ->• oo with the help of estimate (5.6) for |/ m ,fc|. 
We set Q := -R 2 /<? < l/R and use the sequence f3 k := Q fc |/ mj fc|- Using the elementary 
estimate 

IK^fcllp^lKtolli-IK^fclli^lK^fcllfi, 

we obtain 

fe<m (5.13) 

<C||(a fc ) fc ||,2 p- 2m (|^ | + (i? 2 /?)) m - 

To guarantee fast decay of IaC 1 *], we choose e and £i such that 

1 f R 2 \ 1 . . 

e + — <-5- (5.14) 



(1 - eO 2 V 9/ -R 

This is possible since q > R 3 . Combined with our assumptions \z \ < e , and 1 — p < £i, 
we have obtained the estimate 

|/C ut | <C||K) fe || P iT™. (5.15) 

We next study the other error contribution F in in (5.10). Our goal is to express, 
analogous to (5.11), 

F in = AAufclin = ^»W-(^- *o) m ) , (5.16) 

fcGN m>l 

and to provide an estimate for the coefficients /x™. 

As a first step we expand the function Vk(z) = $l(z~ k ) on <9E, which was done in 
(5.5). Since both components of the function z > (z — z ) m are harmonic in E, the 
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expansion of the boundary values provides us also with the harmonic extension v k \x. 
Formula (5.5) yields 

v k ( Z ) = ^Y.^p 2rn+1 y lJ ^ ■ (* - z ^ m forzes ' 

m>k 

We next evaluate the normal derivative, using again (5.1). We find, on <9E, 

m>k P m>k ^ 

This provides for the normal derivative from inside £ the expansion (5.16) with the 
coefficients 

^ = E A ^ 2m+1 r lj ^- • ( 5 - 17 ) 

k<m P 

This expression for u 1 ^ is analogous to (5.12) for In particular, /z™ can be treated 

similarly to in (5.13). 

To sum up, we have obtained 

oo 

F-K^/i^^-^r^ 1 ^, with |/i m | <C||(a fc )klli 2 R~ m - (5-18) 

m=l 



Step lc. Correcting the error on <9£. We now correct the error term F given by (5.18). 
We recall that the error was introduced by V v through V • (AVV V ) = F. 

We define m* to be the smallest integer with (p/ R) 2m * < rj and decompose accord- 
ingly 

F = F low H l L<9£ + F high H 1 [dZ 

m<m* m>m* 

We will correct the high frequency error F hlgh by taking w v to be the solution to Aw, = 

The low frequency error F low must be treated with a quite different approach. The 
basic idea is to use the perfect plasmon waves V k and V k as in (3.3). We define V k (z) = 
*k{z k ) and V k (z) = S(z k ) for z E B R (0), and V k (z) = R 2k U{z~ k ) and V k (z) = R 2k Q( Z - k ) 
for z G M 2 \ Br(0). These functions are perfect plasmon waves for the curve dB R (0), 
but they are not solutions on <9£. In other words, the nonzero functions V • (AW k ) and 
V • (AW k ) are concentrated on <9£. 

The normal derivatives on <9£ of these functions have been used before. There holds 
v ' VVfe — p~ 1 k^R,((z — z^){z k ~ x )\ compare (5.1). Similarly, we have for the imaginary 
part the normal derivative v ■ W k = $l(z — Zq ■ ikp~ x z k ~ x } = p~ 1 A;$5((z — zo){z k ^ 1 )). 

The fact that the functions V k and V k are not solutions on <9X can be used to our 
advantage. We expand the low frequency error _p low in terms of the residuals of these 
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functions. Expanding with respect to the center z = 0, we find, on <9E 

m* m* m , \ 

f 1 ™ = % j2 »m{z - z r-\z -^) = »E^Eu-i) ** -1 (-*b) ra ~*(* - *b) 

m=l m=l fc=l ^ ' 

= E (* X>- (T- 1) (-^) m " fc ) ^ fe_1 (* - ^)) 

k=l \ m=fc ^ ' / 

"E (Til 1 ) (-^) m - fe ) 3(^-^o)) 

fe=l \ m=fc ^ ' ) 

m* m* 



k=l k=l 

with the real coefficients Pk and (3k given by 



TO* / 

h-ih = % E^ uli 



We use now the estimate (5.18), \/i m \ < CR m , to estimate the complex coefficient 

Pk = Pk + 0k: 

m* , % m* 

1/5*1 < c w u Ii ) w™ - * ^ c E i2_m ( 1 + i^i) m_1 • 

We choose e such that 

(l + e )<R. (5.19) 

Hence, as (1 + \z \)/R < 1, we have 

m* 

< C R ~ m ( l + l^o I 1 < CR-\l + l^ol)^ 1 • (5.20) 



m=k 



We can therefore compensate the low frequency errors introduced by V v of (5.9) using 
the functions Vk and V k - We recall that V • AVVk = —2d u Vk'H 1 [dT, and set therefore 

j m* j m* 

^ := n + ^ E ^ + 2 E • ( 5 - 21 ) 

k=l k=l 

With this choice of v v , we have V • (AVv v ) — F + f — F low = f + F high . If we choose w v 
as the solution of Aw^ = F hlgh , we obtain 

V • (AVv v ) - Aw v = f + F high - F high = / . 

In particular, the constraint of (2.10) is satisfied. 
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Step 2. Calculation of energies. It remains to calculate the energy I^v^^w^). By the 
triangle inequality, the energy of v n is bounded if we can control the energy of each term 
on the right hand side of (5.21). 

Recall that Vk defined in (5.8) agrees with v k defined in (4.6) on R 2 \E. Accordingly 
the energy contribution of V v is bounded by a similar argument as in the proof of Propo- 
sition 4.1. But since we do not have orthogonality of the basis functions, we calculate 
here with the ^-assumption on at- Using (5.9), the triangle inequality, and the fact that 
we are in the case q > R 2 , we find 



1/2 / r \ 1 / 2 

(5.22) 



V J |W„| 2 ) = v^||Wj L2 < l Afc l (/ l V{}fc | 2 ) 



k k 



The calculations for the energies related to the corrections involving Vk and Vk are 
identical, we therefore treat here only the contribution of ^2$kVk- Exploiting orthogo- 
nality and the estimate (5.20) for /3 k , we find 



k=l 



2 m* 

<c v f2\Pk\ 2 J \vv k \ 2 



k=l 



< CriJ2^ 2k ( l + \zo\) 2h kR 2k < C(p/R) 2m *(m*) 2 (l + \z \ 



,2m* 



k=l 



where in the last inequality, we have used the fact that r\ < C(p/R) 2m * by our choice of 
m*. By the assumption p < 1 and the choice of Eq in (5.19), the energy is bounded. 

It remains to estimate the energy contribution of w v , given by the solution to Aw v = 
F hlgh . Note that the squared norm of _F hlgh can be estimated by (5.18) as 

HF^H^ < C J2 \»™\ 2 p 2m < C\\(a k )k\\l (p/R) 2m * < Crj. 

m>m* 

By the properties of the solution operator (— A)^ 1 acting on functions in H^ 1 (R 2 ) with 
vanishing average, we conclude that the energy contribution n^ 1 f |Vw^| 2 < Cr]^ 1 n < C 
is bounded. 

The proof of non- resonance is complete by inequality (2.13). □ 

Regarding the last proof we remark that the decomposition of F into high and low 
frequency parts was only necessary in order to obtain an improved lower bound for q. 
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A Are there perfect plasmon waves in dimensions 

The perfect plasmon waves ipk, defined in (3.3), play a central role in the construction 
of trial functions in our variational arguments. In this section we consider the question 
of whether such plasmons exist in dimensions n^2. Plasmon solutions may be viewed 
as solutions to the plasmonic eigenvalue problem (see [7] and references cited therein). 
Specifically, for a fixed radius, R > 0, we seek ip(x), x G M n , such that 

v • (Avv>) = °; i> -»■ as r = \A -> 00 (A- 1 ) 

where 

f 1 \x\ > R- 

A(x) = { ' l x l >/1 ' (A.2) 



This is the situation where we take the core £ = 0. 

A function ip is a weak solution of (A.1)-(A.2) if and only if 

Aip = 0, \x\ ^ R (A.3) 

^\\ x \=r- = ^\\x\=r+ (continuity) (A.4) 



dr 



\x\=R~ 9r 



= (flux continuity) (A. 5) 

\x\=R+ 



ip — > as \x\ — > 00 . (A. 6) 

We are interested in the case e — — 1 and we show 

Claim 1: There are no plasmons localized on spheres in dimension n — 3. 

Claim 2: There are plasmons localized at the plane separating half-spaces in any spatial 
dimension n > 2. That is, for any n > 2, there are solutions of (A.3)-(A.6) in the case 
where the spherical interface \x\ = R is replaced by the planar interface x n = 0. 

Proof of Claim 1: Let Y l (Q), il £ S 2 denote any spherical harmonic of order I. That is, 

-A S 2Y l {Q) = l(l + l)Y l {Q), l>0. (A.7) 

The corresponding solutions of Laplace's equation in dimension 3 are: 

r l Y l (n) and r^ 1 Y l (Q), r = \x\. (A.8) 

Regularity away from the interface and decay at infinity imply 



il>(x) 



r l Y l {Q) , < r < R; 

cr -i-^Y l (Q), r>R. 



The constant c is determined by the interface conditions. Continuity implies c = R 2l+1 . 
Furthermore, the left hand side of (A. 5) is equal to — 7^ 0. Hence there are no perfect 
spherical plasmon waves in dimension 3. A similar calculation holds in all dimensions 
n > 4, where -A S n-iY l (n) =l(l + n- 2)Y l (tt). 
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Remark A.l. Note that 



r l Y l (Q) , 0<r<R 
R 2i+i r -i-i y i(ty , r > R 



is a plasmonic eigenstate with corresponding eigenvalue E\ = —j^i- This sequence of 
plasmonic eigenvalues approaches — 1 as I — > oo. 

Proof of Claim 2: For any n > 2, we write x G W 1 as x = (x±,x n ) and define: 



Then ip(x; £) is a plasmon wave. 
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x n > 
x n < 



where C„eK and £ ± e K' 



are chosen so that: 



\u 2 = a -a- 



